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Abstract 

We investigate brane inflation driven by two stacks of mobile branes in a throat. The stack closest to the 
bottom of the throat annihilates first with antibranes, resulting in particle production and a change of the 
equation of state parameter w. We calculate analytically some observable signatures of the collision; related 
decays are common in multi-field inflation, providing the motivation for this case study. The discontinuity in 
w enters the matching conditions relating perturbations in the remaining degree of freedom before and after 
the collision, affecting the power-spectrum of curvature perturbations. We find an oscillatory modulation of 
the power-spectrum for scales within the horizon at the time of the collision, and a slightly redder spectrum 
on super-horizon scales. We comment on implications for staggered inflation. 
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I. INTRODUCTION 



Constructing models of inflation in string theory is an active field, for reviews see e.g. 
Because string theory naturally has many dynamical degrees of freedom, multi-field models of 



inflation are also common. 
1^12] or multiple D-branes 



ixamples are inflation from axions [7,], tachyons 



M5-branes 



13l4l5l| to name a few. If many degrees of freedom are present, it is 



natural for them to decay or stabilize during inflation one after the other. A phenomenological 
approach to recover some of the effects caused by decaying flelds was proposed in 16Hl8t| and 
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dubbed staggered inflation. [I6l4l8l| relies on coarse graining, assuming that many fields decay in 
any given Hubble time, which may or may not be the case in concrete models. 

In this paper, we set out to investigate a single decay in detail to extract analytically the 
consequences for cosmological perturbations. 



20|427| 



As a simple and yet non-trivial framework, we choose the KKLMMT [3] brane inflation 
set up. In this proposal, inflation is driven by a brane-antibrane pair, located in a warped throat 
generated by the background branes and fluxes 28l43Cll|. Inflation ends when the distance between 



the mobile brane and the antibrane reaches the string scale and a tachyon is formed. After the 
brane-antibrane collision, the energy stored in their tensions is released into light closed string 
modes, which can ignite reheating. We extend this proposal by including two mobile stacks with 
pi and p2 branes respectively and pi +P2 antibranes at the tip of the throat. We assume that the 
stacks are separated from each other so the stack closer to the tip, say stack one, annihilates first 
during inflation. This collision results in particle production during inflation and a sudden change 
in the inflationary potential. We would like to examine some effects of the produced particles via 
the change in the equation of state parameter on the power-spectrum of cosmological perturba- 
tions. Inflation ends when the remaining stack of p2 branes annihilates with the background's p2 
antibranes. This extension does not improve upon the back-reaction and fine-tuning issues of the 
KKLMMT proposal 
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35l |. but provides a well motivated toy model in string theory. 



To handle the analysis, we assume that perturbations in the effective field describing the first 
stack, 501, are taken over by the radiation bath, decaying quickly after the annihilation event. 
Hence, we focus on perturbations in the remaining field only, 5(/)2. For this simplification to 
be consistent, we assume that pi P2- We derive asymptotic solutions in the two slow-roll 
regimes before and after the collision for this field, which we match at the collision time. Treating 
the mechanism of tachyon formation and the collision as instantaneous events, we use the Israel 



junction conditions 



361 ] to provide the matching conditions, which boil down to the continuity of 



the Bardeen potential <I> and the curvature perturbation IZ in the case at hand. It is in these 
matching conditions that the jump in the equation of state parameter enters. We then compute 
the power-spectrum of curvature perturbations and find a slightly redder spectrum on scales that 
were super-horizon at the time of the first stack's annihilation. On sub-horizon scales, we find an 
oscillatory modulation of the power spectrum. 

The concrete outline of this paper is as follows: we review and extend the KKLMMT setup in 
Sec. |TI] before deriving the background, slow-roll solution in Sec. IIIII Perturbations are discussed 
in Sec. \TV\ with the matching conditions in Sec. IIVDI and the Bogoliubov coefficients in Sec. lIVFi 
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These coefficients enable the computation of the power-spectrum in Sec. |Vl which we expand for 
super- and sub-horizon scales in Sec. IV Al and Sec. IV B[ We conclude in Sec. IVII 



II. BRANE INFLATION WITH TWO STACKS 

Here we present our set up which is an extension of the KKLMMT proposal of warped brane 
inflation jiol . 

We start with a warped throat given by the geometry 

ds^ = h-^/^ir)dx''dx'' + h^/^{r){dr^ + r^dsj) (1) 
where h{r) is the warp factor 

MO = ^ . (2) 

Here L is the Ads scale of the throat, which is created by N coincident background branes located 
at r = 

= CQsNa'^ , (3) 



where gs is the perturbative string coupling, y a' is the string theory length scale and c is a 
geometric factor depending on the internal geometry, ^,35. For example, for we have c = Air 
while for T^^'^^ we get c = 27tt/A. Alternatively, one may imagine that the background ([1]) is 



created by turning on fluxes 
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29[. As usual in models of brane inflation, it is assumed that the 
throat is smoothly glued to the bulk of the compactification. 

The metric in ([T|) represents the background geometry, on top of which we add two stacks of pi 
and p2 D3-branes within the throat at positions ri and r2 respectively. At the bottom of the throat 
are pi + p2 D3-branes. In accordance to the KKLMMT proposal, the inflationary period arises 
from the mutual attraction of the mobile D3 branes and anti D3-branes plus the interaction of 
the D3-branes with the background fluxes. We work in the probe brane approximation, meaning, 
we assume that the added stacks of branes and antibranes do not destroy the background. For a 
schematic view of our set up see Fig. [TJ 

A crucial feature of our model is the annihilation of brane-antibrane pairs during inflation. To 
illustrate, during inflation, the mobile stack of branes closest to the tip of the throat, comprised of 
pi D3-branes, annihilates with pi D3-branes. The effective field corresponding to this stack drops 
out of the dynamics, leading to a sudden jump in the potential. After this annihilation, inflation 
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P2 D j-branes 
Pi D3-branes 

To (pj + Pj) Dj-branes 

FIG. 1: A schematic view of the set up. Two stacks of mobile D3-branes are located at ri and r2 while 
D3-branes are located at the bottom of the throat at tq. The first annihilation happens when the physical 
distance between ri and ro becomes equal to string scale Is, which occurs shortly after rji — —1. 



proceeds driven by the remaining stack of branes. Inflation ends when the second stack annihilates 
with the p2 antibranes located at the tip. It is straightforward to generalize this setup to multiple 
stacks, Pi, I > 2, and multiple collisions. However, since our primary goal is the investigation of a 
single collision event, we consider only two stacks. For a schematic of the collision, see Fig. [2j 

After the annihilation event, the kinetic energies of the branes as well as potential energies 
originating from the branes tensions are transferred into closed string modes ^ . In order to simplify 
the analysis, we assume that these are in the form of massless modes so that they behave as 
radiation. As we shall see in our analysis, it will not change the results if one assumes that the 
produced particles are in the form of dust. The important effect is the sudden change in the equation 
of state parameter due to the collision. Furthermore, we assume that the brane annihilation and 
tachyon formation is sudden, leading to an efficient method of converting decay products of the 
brane-antibrane annihilation into closed strin g m odes. For a list of phenomenological studies of 
reheating in models of brane inflation see e.g. 
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471 ] and the reference therein. 



^ We ignore other decay channels in this study, such as the productions of cosmic strings [4^ . . Our conclusions 
primarily depend on a jump of the equation of state parameter; the magnitude of this jump would not change 
much if other, subdominant decay products different from closed string modes were included. In their presence, 
we expect only small quantitative changes in i.e. the power-spectrum after the collision. 
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Annihilation Aftermath 
t=tc 




FIG. 2: A schematic view of the annihilation event at tc- Once the pi D3-branes come within Ig to the 
antibranes, they annihilate. Their decay products, such as light closed string loops denoted by pr, cause a 
jump in the equation of state parameter w (see Sec. lIVEp . We require p2 S> pi in this study, so that is 
subdominant after the collision. 



To calculate the inflationary potential, we treat the pi +p2D3-branes as probes in the geometry 
created by the stacks of N background branes as well as the additional two stacks of D3-branes 

m 



3l\ . This leads to the condition pi^2 ^ N. The warp factor in this configuration is then given 



by 



Kn = — + 



+ 



P2 



(4) 



r"* A^|r — ri|^ N \r — r2\^ 
This form of the metric indicates that the multiple stacks of D3-branes in the absence of antibranes 
and fluxes are supersymmetric. 

The action of the pi +p2D3-branes located at the bottom of the throat, rg, is [lol 

S = -2(pi+p2)T3 I d^xh{rar\ (5) 

where T3 = l/{2'K)^gga''^ is the D3-brane tension and the factor of 2 in Eq. ([5]) originates from the 
combined tension and charge of the antibranes in this background. 



Defining = \/~piTlfi for / = 1^ (pA = \l (pi + P2)?3''o where the subscript A denotes 

3, 



antibranes, and noting that ro <C ri^2 
before the collision as 



we can write the potential for the two scalar fields 4>i 



V- 



where 



2{P1+P2)T3^ 



p](j)[ 



N{pi +P2f 



(6) 



(7) 



If ri < r2 the field (pi drops out of the dynamics first. After the annihilation of the first stack, 
the potential reduces to 



where 



1 - — 
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(8) 



^0^ = 2p2r3-^, (9) 

leading to a jump in the potential. 

After coupling the system to four-dimensional gravity, the total action is 

1 fl \ 

-^R - -9^(/.2a^(/'2 - ( i^d^,(t>id^'(t)i + v-{4>u(t>2) J e{-t + tc) 



S = J d x\J —g 

+ y d^x^g {Cra - V^+(02)) e{t - Q (10) 

where Mp is the Planck mass related to Newton's constant G via Mp = l/(87rG) and tc indicates 
the collision time. In order to respect conservation of energy at tc we are forced to include an addi- 
tional contribution to the action, Cra, corresponding to radiation created during the annihilation 
event. 

For self-consistency, we require the initial distance between branes and antibranes to be greater 
than the string scale. Is = rnj^, otherwise, the appearance of tachyons leads to the immediate 
annihilation of branes and antibranes. The first inflationary stage terminates at t = tc where the 
distance between the first stack at ri and the antibranes at the tip is Is- Defining (f)i{t = tc) = 4>ic 
and using the metric in ([T]), we obtain 4>ic = (/'a exp(/s/L). The onset of tachyon formation 
and brane annihilation results in the violation of the slow-roll conditions at the collision time. 
Practically, as in [49j|, the onset of the brane collision is indicated by f/]"(tc) = —1 where r/i, defined 
in p4p . is a slow-roll parameter. 

III. BACKGROUND: SLOW-ROLL INFLATION 

In this section, we derive the background dynamics in a homogeneous and isotropic universe 
where the four dimensional metric is 

ds2 ^ _^^2 ^ a{tfd^2 ^ ^^^^ 

a{t) is the scale factor and t is cosmic time. The Friedmann and Klein Gordon equations are 

H^=P, (12) 
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dV 

cl)j + 3H^i + — = 0, (13) 
d(pi 

where p is the total energy density containing contributions of scalar fields and radiation, pr, after 
the collision. 

The background solutions for before the collision are easily obtained if the system is in the 
slow-roll regime, that is if the slow-roll parameters 

= , m^Mj,^ , e^e, + e, , e,, = ^2 , (14) 

are all much smaller than one. Since there is no cross term in the potential, the slow-roll parameter 
r/12 oc V0i(^2 vanishes. Here, we used the short-hand notation = dV/d4)i. During slow-roll we 
can approximate 

3//fe,(^/ + ^^^0, (15) 
where we used the potential in ([6]) and 



is the expansion rate in the slow-roll limit before the collision. Eq. (jlSp can be solved to 

^f.n-^fc-%^*e, (17) 
^bc 

where we set the initial time at the start of inflation to zero and </>/(t = 0) = 4>iin as well as 
4>i{t = tc) = <pic- We can replace cosmic time in terms of the number of e-folds dNf. = Hdt, so 
that 

</.f - </.f , 24M|, bi (Nt - AN,) . (18) 

Here Nt — 60 represents the total number of e-folds before the end of inflation, which needs to be 
about sixty in order to solve the flatness and horizon problem, and AA'^g is the number of e-folds 
from the time of the first stacks's collision till the end of inflation, which we assume to be close to 
Nt throughout this article. After the first stack's collision, the field 0i drops out of the dynamics 
and the system consists of the field (/)2 and radiation, which is diluted quickly. Thereafter, (f)2 
resumes its slow-roll evolution, that is the dynamics of (p2 is given by 

3Haj2 + ^^^0, (19) 
Y2 
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where Hac is the expansion rate in the slow-roll limit after the collision, 



Hac^yj^. (20) 



Eq. p9j) is solved by 



'Ac-<i^2f- 24M|, 62 AA^e • (21) 



Here (^2/ is the final value of (/)2 when inflation ends, corresponding to the time when the second 
stack annihilates with the remaining p2 antibranes at the bottom of the throat. 

The first stack's collision is instigated when i]^{tc) = —1, as in j^], where rjf represents the 
slow-roll parameters before and after the collision. The rjj can be calculated to 



rji = -20- 

cpj 

(22) 



20M2 (ro 



NT^rl \rj 

Consequently, the condition ryf (tc) = — 1 results in 

(j)\^ = 2GbiMl , (/.?i„ = 246iM|,(iVT-A7Ve + 5/6). (23) 

Similarly, tj, the time of the second stack's collision and the end of inflation, is set by 'i]2 i^f) ~ 
which results in (j)^ ~ 20b2Mp. Combining this with (j2ip and (jlSp we obtain 

06.^~24M|,62(iVT + 5/6) , (/.^,~24Af|,62(AiVe + 5/6). (24) 

For later reference it is instructive to express the slow-roll parameters in terms of the number 
of e-folds AA^e and Nt- Using Eqs. ([23]) and we get 

Here ry^-^ and 77^-^ indicate the slow-roll parameters at the start of inflation whereas ry^^ indicates 
the corresponding slow-roll parameter after the first stack's collision. To obtain 77^^ we assume 
that radiation is quickly diluted so that the remaining stack resumes its slow-roll motion quickly. 
Regarding the other slow-roll parameters, we note that e/, e <C since 



e/ 2pj ( ro 



|r//| 5A^ \ri 
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- «1, (26) 



due to t^ 
such as 



le large value of N ^ pj and < rj. This is common in warped brane inflation models 
I9I . since the warp factor reduces the height of the potential. This implies that £/ <C |r?j| 
which is the main reason why gravitational wave production in these models is suppressed 13]. 



Since (pi drops out of the dynamics after the colhsion, the energy associated with the tensions of 
the pi antibranes as well as the kinetic energy of the pi mobile branes are transferred to radiation. 
During the annihilation event the velocity 02 is continuous, but due to the step in the potential 
([8]), the acceleration (j)2 makes a jump. Energy conservation, V^T^*^ = 0, dictates the value of pr 
just after the brane annihilation, 



Pr{tc 



(27) 



where we defined [f]± = f{t'^) — fit^)- Assuming slow-roll inflation before the collision, one can 
neglect the kinetic energy of the mobile branes compared to the tension of the antibranes and 
obtains 

Pr{t,) ~ {v^ - v^) = 2piTs 4/1" . (28) 

We are interested in the limit where the energy transferred into radiation is small compared to the 
background inflationary potential, that is Pr{tc)/vo = Pi/{pi + P2) ^ 1 or pi ^ p2- 

Next, we would like to solve for the evolution of the scale factor after the collision; noting that 
Pr scales like radiation 



Prit) 



(29) 



° -^^Va(t), 

where Oc is the value of the scale factor at tc, we can approximate the Friedmann equation for 
t > tchy 

/ ■ \ 2 

,+ ^ /„- _ .,+\ I \ (30) 



which can be solved analytically for a(t), 

1/2 



a{t) 



1 + \ Vq /vq 



1 



1/2 



Hac{t — ic) 



(31) 



(1 + ^f^r 

As explained before, Hac, given by Eq. (j20p . is the expansion rate after the collision, when radia- 
tion is diluted and the approximate de-Sitter background is recovered from the interaction of the 
remaining p2 D3 and D3-branes. Comparing Hac with the expansion rate before the collision, we 

Hhc 7^ Hac due to the presence of radiation although 



(32) 



get Hac ^ Hi,c\Jv^ /vq . Note that H{t+' 

H itself is continuous. Asymptotically, the scale factor in ([3T|) approaches 

1/2 



a{t > tc) 



1 + V /vq 



ace 



An interesting implication of (j32p is the increase of the scale factor's amplitude by a factor of 



{1 + y Vq /vq)^/"^ /V^, caused by the transient presence of radiation. 
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IV. PERTURBATIONS 



In this section, we consider linear perturbations before and after the colhsion at t = tc- For 
t < tc perturbations are carried by the two scalar fields (we assume a flat field space metric) and 
after the collision, by the remaining scalar field and radiation. Perturbations in 0i are carried over 
predominantly by the radiation bath for t > and decay rapidly in the subsequent inflationary 
phase. Thus, we do not need to follow their evolution. 

Our goal in this section is to derive and solve the relevant equations of motion for perturbations 
in the surviving field before and after the collision, which need to be matched at tc- 



A. Equations of Motion 



The most general line element including scalar perturbations is 



ds'^ = -(1 + 2A)dt^ + 2aB^idx'dt + a^[(l - 2^/;)5ij + 2E^ij]dx'dx^ 



(33) 



where we did not remove any gauge modes yet. The equations of motion of the 



the scalar fields include couplings to the metric degrees of freedom and read 



51 



perturbations in 



52| 



bj + 3H6(l)j + —5(j)i + V<i>j<t>j6^j = -2V^jA + 0/ 



k' 

A + Sip + ^l a^E - aB 



(34) 



with / = lj_2. One can define gauge invariant metric perturbations, for instance the two Bardeen 



potentials 



5C| 



^ = A + (aB-a'^E) , 



^ = ip - H (aB - a^E) , 



(35) 
(36) 

that coincide in the absence of anisotropic stress <I> = Further, we can define the gauge invariant 
Sasaki-Mukhanov variables 



Qi = ^4>i + 

li 



(37) 



These coincide with the field perturbations in the spatially fiat gauge 1^ = 0. In the absence of 
any other components in the energy momentum tensor, such as radiation, and usin g th e perturbed 
Einstein equations as well as the background equations of motion one gets (see 5l|, [53] for details) 



(38) 
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Thus, the equations of motion for the Sasaki-Mukhanov variables decouple from the perturbed 
Einstein equations for the Bardeen potentials, and we need not be concerned about the latter. In 
deriving (|38p we neglected perturbations in the radiation bath. Our philosophy is to keep the setup 
as simple as possible, while retaining some crucial effects onto perturbations of the annihilation 
event. 

Following [53], one can show that the equations of motion simplify during slow-roll to 

where a prime denotes a derivative with respect to conformal time d/dr = ad/dt, ui = aQi and 
the matrix Mjj is defined as 

[ e + '^^i-mi 2ei2 \ 

Mjj=\ , (40) 

y 2ei2 e + 2e2 - 7722 / 

where we used the slow-roll parameters defined in ()14p ^. 

B. Before the Collision 

Since e/ <C in warped brane inflation, the matrix Mjj becomes diagonal and the equations 
for the Qj separate, 

u'l+l^k'-i^^^^j^O, (41) 

where 

^7 = ^ - ^7 (42) 

and rjj is defined in Eq. ([22]) . 

We impose the Minkowski vacuum state ui — e~^^'^/V^ in the far past (r —00), corre- 
sponding to the choice (see i.e. the review js^ ) 

uj = :E^e^A>^7+im/2H^^}(^_kr) e, , (43) 
2 f-i 

where H^^l{x) is the Hankel function of the first kind of order and ei are independent unit 
Gaussian random fields with 

<ej>=0 , <e7(k)ej(k') >=(5/j53(k-k'). (44) 



^ We use the potential slow-roll parameters whereas [531 uses the Hubble slow-roll parameters, which coincide at 
linear order only. 
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Since perturbations originate during the early stages of inflation well before the collision, we can 
use rji = r]J-^ so that fij = 3/2 — -qj-^ from ([2^ . resulting in 

_ _ 3 1 _ _ 3 1 . 

^1 -2 + l + |(iVT-AiVe) ' ^2-2 + 1 + 6^^- ^""^^ 



C. After the Collision 

After the collision, when radiation is diluted away and (j)2 resumes its slow-roll motion, the 
equation for is given by ([39]) with r/^ — )• 77^ and fi2^f^2^ which can be integrated to 

u+ = ^^e-(^^+i/2)/2 (aH^'}{-kT)+(3H^'}{-kT)) . (46) 
2 V ^2 / 

Here /U^ = | — f/J, while a and /?, often referred to as Bogoliubov coefficients, need to be determined 

by an appropriate matching procedure near the collision. Using Eq. ()25p we have 



= - + i . (47) 

2 1 + fAiVe 

We note that since Nt ~ AA'^e ^ 1- 

In deriving the asymptotic solutions in (I43p and (I46p we treat the slow roll parameters as small 
and constant. This approximation is common when dealing with inflationary slow-roll models and 
its validity has been shown numerically. In our model, slow roll is briefly violated during the 
collision, but due to the short time scale of this collision, we expect this approximation to remain 
valid. The time scale of tachyon formation and brane and anti-brane annihilation at the bottom of 
the throat is given by the inverse of the warped string scale. To trust our four-dimensional effective 
field theory approximation, the Hubble expansion rate has to be much smaller than the warped 
string scale. This indicates that branes annihilate and slow-roll is violated for less than an e-fold 
and the slow-roll approximation is restored quickly. 

Our goal in the next sections is to calculate a and f3 to compute the power-spectrum. For later 
use, we define 

a = aiei + 0:262 , (3 = f3iei + (32^2 ■ (48) 
D. Matching Conditions at tc 

We treat the collision in a sudden approximation, infusing the entire energy of 0i into the 
radiation bath pr at t^- We are interested in the limit where the energy transferred into radiation 
is small compared to the background inflationary potential and pr{tc)/vQ = pi/{pi + P2) ^ 1- 
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To perform the matching of cosmological perturbations, we demand that both the intrinsic and 
extrinsic curvatures are continuous on the hyper-surface separating the two inflationary phases 



36l . l55l . l56l |. As in standard hybrid inflation models [57l . |58 | the hyper-surface of the first stack's 
collision, and thus the transition between the two phases, is set by a certain value of an inflaton 
field, (f)i = (f)c for us. Treating the transition as instantaneous, that is assuming At ^ H^^(tc), 
the cosmological matching conditions relating perturbations before and after the collision become 



m± = , [7^]± = 0, (49) 

where the absence of anisotropic stress was used (<I> = ^'). These matching conditions are valid for 
modes with k <^ (At)^^, which includes modes that are within the Hubble horizon at t^- In the 
above, the comoving curvature perturbation is defined as 

2 \ 

The matching conditions in (|49p for perturbations are accompanied by the continuity of the scale 
factor and the Hubble parameter T-L = a' /a (and thus the total energy) at the background level, 

[a]± = , [^]± = 0, (51) 

which also follow from the continuity of the induced metric as well as the extrinsic curvature on 
the hyper-surface set by the first stack's collision. 

To make contact with observations, we need to compute the power-spectrum of the curvature 



perturbation TZ, which is related to the Sasaki-Mukhanov variables via 

7^ = 



5li 




3(l + u;)M27^^^^^^' ^^^^ 

if no additional contributions to the energy momentum tensor are present. 

To impose the matching conditions, we need to know (<I>, TZ) as a function of the Qj. Fortunately, 
TZ is already given in terms of the Qj in Eq. (j52]) . To obtain $ as a function of the Qj , we use the 
perturbed Einstein equations [50| 



(^^ + ^0. = (54) 



p 
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with the gauge invariant perturbations of the energy momentum tensor 

= ^ ^ (-0'/^ + ) ' + V,^aHct>f^) , (55) 



and the gauge invariant field perturbation 



(t>i + (t)'i{B - E') . (57) 

In the flat gauge we have if) = Q and we can use Qi = Scjii as well as 6(j)\^^^ = Qi — (p'j^ /H. 

Plugging these into the perturbed Einstein equations and using ^' = <^ as well as the background 
equations of motion, we can solve (|54p for <!>', insert it into (|53|) and arrive after some algebra at 



w)(i)'jn 



'iQ'i , (58) 



where we used the slow-roll approximation and worked in leading order of the slow-roll parameters, 
that is we used 1 -|- ~ and a'^V(f,j + ^(p'j'H ~ before and after the collision. 

It is at this point where we make an important approximation to simplify the matching proce- 
dure: based on the fact that perturbations in 0i are predominantly carried over by perturbations 
in pr, which rapidly redshift, and that p2 ^ pi so that pr <^ p^j, we ignore Qi before the collision 
just like we ignored 6pr after the collision. Thus, in order to perform the matching, we approximate 
(f52]) and §^ by 

-^'^^^'^2^2, (59) 

n w — ^(t>2Q2 , ■ (60) 



0, (61) 



so that the matching conditions in (|49p become 

Q2 



[Q2]± = o 



.1 + 

where we used that ^'^ is continuous. If we now use the asymptotic solutions before and after 
the collision respectively, ()43p and (|46p . we can read off the Bogoliubov coefficients a and /3. Our 
approximation of neglecting perturbations in radiation and is equivalent to having a ~ 012^2 
and ^^h^2 m Eq. 

Before we proceed, we would like to review our approximations. Our first assumption is that the 
system is in the slow-roll limit before the collision. A further crucial assumption in our analysis is 

15 



that radiation is quickly diluted after the collision so that the system resumes its slow-roll evolution 
quickly. In this limit, (|5'2p and (|58p can be used, which are valid for systems composed of scalar 
fields only. We note that there have been many phenomenological studies of brane annihilations 



and reheating in models of brane inflation e.g. 



401445(1 . However, we do not know the exact stringy 



details of tachyon formation and the brane collisions, nor the mechanism of energy transfer into 
radiation, especially in a non-trivial background with background fluxes present and collisions 



6^, 



63] • Nevertheless, we expect that 



between multiple branes and anti-braens, see for example 
perturbations in are predominantly carried over by perturbations in radiation so that we can 
use ()59p and ()60p in order to perform the matching. We still keep some effects of radiation via the 
change in the equation of state parameter w; by means of the Bogoliubov coefficients we can relate 
the solutions in the two slow-roll inflationary regimes by using the matching conditions (j6ip . One 
shortcoming of this treatment is that we do not take into account perturbations in radiation and 
01, which would enter into the matching conditions. Nevertheless, even with this simplified setup, 
we find interesting effects that we expect to prevail in a more rigorous treatment, which we leave 
to future study. 



E. The Equation of State Parameter 

Before the collision, the branes move slowly and the slow-roll conditions are satisfied for both 
fields so that W- ~ —1. Close to the annihilation event, |r/j~| becomes of order one; however, since 
the annihilation takes place shortly thereafter, the field (/>i has no time to speed up and we may 
use its slow-roll value, even though slow-roll is violated. Particularly, we can neglect its kinetic 
energy compared to the potential one just before and after the collision. Similarly, we may use the 
slow-roll value for 02- 

We would like to relate the outgoing solution in the slow-roll approximation in ()46p to the 
incoming solution in ()43p . We already specified the parameters ^'J and /i^ in (j45p and (j47p . but 
we still need to calculate the equation of state parameter for the incoming and outgoing solutions. 

For t < tc we have 

where i?f,c is the Hubble expansion rate in the slow-roll approximation from (jl6p . After the 
collision, assuming that radiation is diluted quickly, the slow-roll conditions are satisfied again and 
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the equation of state parameter is given by 



so that 



Here we defined 



1 + U'H 



(64) 



5/6 

02 (Tc) V P2 



where we used (|23p and (I24p . Similarly, the fractional change in u) is 

Aw _ — 7^ 



(66) 



1 + W- 1 + W- 1 + 7^ 
In the limit pi — )• 0, corresponding to no brane collision during inflation, or P2/P1 ^ i^Ne)^^^, we 
obtain 7 — )• 0, — w_ and Au; — )• as expected. 

Finally, we note that during slow-roll inflation T-L ~ — I/t. In our approximation two slow-roll 
inflationary regimes are glued to each other with a continuous Hubble expansion rate so there is 
no jump in conformal time thus r is continuous on the hyper-surface of the first stack's collision. 



F. The Bogoliubov CoefRcients 

At last we are in a position to match the two phases of inflation and calculate the Bogoliubov 
coefficients a and f3, which we need in order to compute the power-spectrum after the flrst stack's 
collision. Starting with the matching conditions (I6ip we obtain 



1 + 



Q2I 
-Q2I 



Q2I 
Q2I 



(67) 
(68) 



1 + 

In the limit where perturbations in radiation and dcpi are ignored, we have a ~ 02 62 and /3 ~ /32 62 
in Eq. (gHD. With Qj = ui/awe flnd 



Q2 



132 



1 + W- H 1^2 V H 1^2 



1 + W- ^^2 1^2 



(69) 



^ To be precise, t~ /t^ ~ 1+37^(1+to_)/3, so that the difference between r~ and is not only slow- roll suppressed, 
but also suppressed by 7 ; as we shall see shortly, 7 <g f is needed to guarantee that additional features do not 
dominate the power-spectrum. 



17 



where we defined 



X = —kTc , (70) 
^ = |(^2 - Z^^) = 1(^2" )• (71) 



To derive (|69j) we used the identity 

<^Hi'^-HiPH'W = -^. (72) 

^Prom here on, all Hankel functions depend on x and a prime denotes a derivative with respect to x. 
Since kc = — 1/tc is the wave number for which perturbations leave the Hubble horizon at the time 
of the collision, we see that x = k/k^ smaller or bigger than one discriminates between super- and 
sub-horizon modes. With the above coefficients and after replacing the Hankel functions in terms 
of Bessel functions as well as using the corresponding identity for Bessel-functions, J'^{x)Y^{x) — 
J^{x)Y^{x) = —l/iTx, we obtain 



(a*(k')a(k)>-(/3*(k')/3(k)> (|a2p - |/32p) '^^'(k - k') (73) 

'-6^{k-k') (74) 



1 + w. 

-T^^'(k-k'). (75) 

We note that |ap — / 1, since the perturbations in 0i encoded in ai and /3i are "lost" at 
the collision (we assume that they are taken over by radiation and redshifted away rapidly). In 
the limit 7 — t- we recover the usual normalization of the Bogoliubov coefficients. However, for 
arbitrary values of 7 we checked that by including the contributions of ai and f3i the relation 
|a|2 _ |/3|2 = 1 holds. 

V. THE POWER SPECTRUM 

We are interested in the power spectrum Vtz of the curvature perturbation TZ, which is defined 
by (see 



54( 1 for a review) 



3 



S^k - k')Vn = j^< 7^(k')*7^(k) > (76) 
To relate TZ to the outgoing perturbations after the collision in ()46p we use (|52p . yielding 

n ^ S-Q2 . (77) 

92 
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At late times {—kr — )• 0) we can neglect the decaying part of the Hankel functions in (j46p . that is 



P2 



M2 



vr 



kT\'^^ 



2 J 



Consequently, we have 



which results in 



vr V 2 



-M2 



(/32 



a2j 



V27r<^2r(3/2)y 



-kr 



1/32 



02! 



(78) 



(79) 



(80) 



The overall amplitude Vti — {H"^ /2-K(j)2)'^ around sixty e-folds before the end of inflation is set by 
the COBE normalization V-ji ~ 2 x 10~^, which in turn sets the Hubble scale during inflation. 

The important feature in this expression is the transfer function |/32 — a2|^i which encodes some 
effects of the brane annihilation in terms of the Bogoliubov coefficients. The scalar spectral index, 
defined as — 1 = d In "P?^/ din A;, becomes 



1 = 712 - 1 + 



d In 1 02 — /32 1 



where n2 



din A; 

4 — 2/i^ is the usual slow-roll result given by 

5 



n-2 



1 



(81) 



(82) 



As in [19] and [49] one can tune the parameters of the model (A^, g^, rg/i? and mg/Mp) to guarantee 
a sufficiently long inflationary phase and n2 — 0.98. 

Compared to standard single field inflationary models, the power-spectrum is altered by the 
transfer function |a2 — /32| ■, which encodes some of the effects of the brane collision and particle 
production during inflation. Additional contributions to the power-spectrum and higher order 
correlation functions originate via back-scattering of the produced particles onto the inflaton con- 



651 ]: to estimate this effect, which can dominate the power- 



densate, leading to IR-cascading 
spectrum [64] and is not retained in our approach, a better understanding of the brane annihilation 
process is needed. 

With 02 and [^2 from (f69l) . we get 

1^2 



97r^ 



|a2 - ^2! = -TT 



16 Vl + i^- 



M2 



J2+ ( j2_ + y2 



M2 



(83) 



+- 



4 l + w- 
2 



J^+{J-J'-+Y-Y'^] 

tlj ^ fJ-2 M2 ^2 tJ'2 



1 +^^+ J j> 



j2_ + y2 



^^2 



M2 



X 



TT 



J„+J'- 

M2 f^2 



l + w- ^2 



l+W- t^t '"2 
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FIG. 3: The transfer funetion |a2 - /32p from (|H3) plotted over X — k/kc for N^. — 60, ANx — 58, pi — 1 
and P2 = 2.5 x IQ-*. 



We are interested in the shape of the transfer function in different Umits, namely x ^ 1, x ~ 1 and 
X ^ 1, corresponding to modes leaving the horizon before, during and after the collision. The full 
expression is plotted in Fig. [3] , but it is instructive to derive approximate analytic expressions of 
the transfer function in the limiting cases x <^ 1 and x ^ 1. 



A. Super-Hubble Scales, a; ^ 1 



The case x <C 1 represents perturbations which cross the Hubble radius long before the collision, 
k <C kc- In the absence of entropy perturbations these modes are constant during the two slow-roll 
phases before and after the collision and we do not expect strong changes as a result of the collision. 
To see this, we use the small x limit of the Bessel functions in (183p . which results in 



|a2 - /?2| 



r2(^+ + i) V2, 

r2(/i-) /x\2(/^J-M2-) 



3Aw 



r2(^+ + i) ^2 



4(1 + w- 
4(1 + 72) 2 



- , +1 + "^+ 
2V^2+M2YT^ 



1 2 



l + Y 



(84) 
(85) 
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With from ([15]) and from ([T7|) we get 

- (1 + f AATe) (1 + INt) ~ 6 ATtAAT, « ' ^'^^ 

3 1 

_ 2 + l+6jVT/5 _ _ 5 iVT - AA^e _ 

where A'^r ~ 60, 1 ^ AA/'e < A't- is the number of e- folds from the colhsion until the end of inflation 
and we expanded in the last step, see Fig. |3] for a plot. We note a mild suppression on super- 
horizon scales proportional to x'^^^2~^^2 ). This is in agreement with the finding of 



661, 



671], whose 



authors investigated the consequences of a discontinuity in the second derivative of the potential 
in a single inflaton model, which also caused a jump in /i. The scalar spectral index in ()8ip picks 
up an extra contribution on super-horizon scales, 

/ + -\ 5Nt- ANe 

n, - 712 + 2 (/i+ - ^2 ) - ^2 + g ^^AjYe ' ^ ^ 

which is small compared to the leading order slow-roll contribution. In the absence of a collision, 
that is for = ^2 and 7 = 0, we recover \a2 — /32|^ = 1, as expected. 
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B. Sub-Hubble Scales, x > 1 



In the limit x S> 1 we obtain 



^2 ^ 



1 - 27^ sin^ I X 



1] + ( ^ + ^^ 
4 / V 1 + 



Svr 



sm 



/^2 ^ _ ^ 
2 4 



(89) 
(90) 



where we expanded for 7^1. Again, we recover |a — /3p — t- 1 in the Umit 7 — )• 0. We observe that 
within the horizon, pertmbations are modulated and the power-spectrum picks up an unsuppressed 
oscillatory mode with amplitude 



27^ ^ 2^ 

P2 



6 

1 + -AATe 
5 



(91) 



see Fig. O If AN^, ~ ~ 60, these oscillations are within the observational window of CMB 
experiments. In order to prevent them from dominating the power-spectrum, we require 27^ < 0.1. 
As a consequence, we need P2/P1 > 2.5 x 10^. On the other hand, the mobile stacks of pi and p2 
branes are assumed to be probe branes. For consistency, we require large enough background branes 
(flux charges) corresponding to > 10^. This can be easily obtained in the light of developments 
in string theory flux compactifications [3], although, the large hierarchy p2 > lO'^pi may look 
problematic. However, if several branes annihilate successively, we expect these oscillations to be 
averaged out to some extent, alleviating this bound. Thus, it would be interesting to extend the 
model examined in this paper to include several stacks of branes that annihilate one after the other. 
Such a setup falls within the framework of staggered inflation, as proposed in [igI. fl?!. We plan to 



investigate such a model numerically in a forthcoming publication 



to bridge the gap between 



the analytic results of this paper, valid for a single annihilation event, and the analytic results of 
Q) Q]) which are expected to hold in the presence of many annihilation events. 

Small oscillations on top of the power-spectrum, as in Fig. \5\ have recently caught the attention 



of the community 



to 



67 



69l472t|. since it is possible to mimic outliers in the CMB power-spectrum 



73h75I1 . In 1661 . |67I| the ringing pattern is damped with increasing x whereas in our model, similar 



71 



72l ]. we find no such damping. This indicates that all modes within the horizon are 
affected by the annihilation event, not only those with close to the horizon scale. However, 
on the physical grounds one may expect that very small scale modes which effectively live in the 
Minkowski background should not feel the effect of brane annihilation and particle creation. Our 
result of Eq. (j89l) with constant amplitude power spectrum modulation may be due to the sudden 



brane annihilation approximation used in our analysis. In practice, the onset of tachyon formation 



22 



1.02 




large x approximation 



U.BHH 1 , 1 , 1 , 1 , 1 

2 4 6 8 10 

X 

FIG. 5: The transfer function \a2 — /32p in (|83p . solid line, is compared to the large x limit in dashed 
line, for = 60, AiVr = 58, pi = 1 and p2 = 2.5 x 10^. 

and brane annihilation happens at a time scale given by the inverse of the warped string scale. 
This time scale is short enough, i.e. much shorter than H^^, so our approximation for super- 
horizon perturbations is valid. However, for wavelengths shorter than the inverse of the warped 
string mass scale the transition should be smoothed out. As a result, this can remove the power 
spectrum modulation for wavelength shorter than the inverse warped string scale. 

A word of caution regarding the scalar spectral index in ()8ip might be in order: if one inserts 
the transfer function (j83p into (jSip and evaluates for large x one finds a linearly increasing mode 
- this mode is an artifact arising from the attempt to describe an oscillatory power-spectrum by 
means of a simple power law, which is inappropriate. 

C. Robustness of Results 

In order to derive the power-spectrum analytically we have made a series of approximations. One 
important assumption in our analysis is the neglection of perturbations in the decaying field and 
radiation, which would both enter in the matching conditions. We plan to relax this approximation 
in a forthcoming publication - preliminary results show the same oscillations in the power-spectrum 
with a comparable amplitude (up to factors of order one), resulting in a similar bound on the 
number of branes. We ignored other decay channels, for instance the production of cosmic strings 
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at the collision. The oscillatory corrections to the power-spectrum that we find are caused by the 
jump in the equation of state parameter and the potential. These jumps are present regardless of 
the exact type of decay products; hence, we expect only minor quantitative changes to the power- 
spectrum if they were included. We also worked in a sudden decay approximation, matching two 
asymptotic solutions at the collision-time. In reality, the decay takes some time which smoothes out 
all discontinuities while slow roll is briefly violated for the decaying field. As explained previously, 
the time scale of tachyon formation and brane annihilation is controlled by the inverse of the 
warped string scale which is much smaller than H"^ during inflation in order for our effective four- 
dimensional approximation to be valid. In conclusion, our result for super-Hubble perturbations 
with X < 1 is trusted. However, for wavelengths shorter than the inverse of the warped string scale 
the transition should be smoothed out. This in turn may remove the unwanted power spectrum 
modulation for very small scales. 



VI. DISCUSSION AND CONCLUSION 



Motivated by the recent interest in inflationary models with more than one dynamical degree of 
freedom, we investigated the consequences for cosmological perturbations of a decaying field. As 
a model we chose one of the best understood models of inflation in string theory, the KKLMMT 
proposal, and extended it slightly: instead of driving inflation by one brane that annihilates with 
an antibrane at the end of inflation, we allowed for two stacks of branes with pi and p2 branes 
respectively. These two stacks correspond to two inflaton fields at the level of four dimensional 
effective field theory. If the two stacks are separated initially, the first stack can annihilate with 
antibranes during inflation, creating closed string loops (radiation) in the process (we ingored other 
decay channels) . We focused on P2 ^ Pi i in which case the produced radiation does not interrupt 
inflation, but causes a slight jump in the equation of state parameter. This setup has all the 



shortcomings of t 
the inflaton field 



re orig inal KKLMMT proposal, i.e. the back reaction of the volume modulus on 
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351 ] , but it provides a well understood framework to investigate the effects of 



a brane annihilation and particle creation in detail. 

At the perturbed level, we focused on fiuctuations of the surviving infiaton field, with the as- 
sumption that perturbations in the field corresponding to the first stack of branes are taken over by 
the produced radiation, which subsequently decay rapidly. This is our main simplifying assumption. 
We also ignored contributions from IR-cascading. The asymptotic solutions for fluctuations in the 
surviving field during the two inflationary regimes, before and after the collision, were matched at 
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the time of the cohision according to the Israel junction conditions. It is in this matching procedure 
that the jump in the equation of state parameter enters. The resulting power-spectrum contains 
a transfer function, determined by the Bogoliubov coefficients. On scales that are super-horizon 
during the collision, we find a slightly redder spectrum. On sub-horizon scales, the power-spectrum 
is modulated by oscillations, or a ringing, with an amplitude of ~ 2.5 x 10^pi/p2- Such a ringing 
pattern might help explain glitches in the CMB power-spectrum. Since these oscillations should 
not yield corrections bigger than a few percent, we deduce that the remaining stack of branes needs 
to be comprised of many more branes than the first stack, p2 ^ 2.5 x lO^pi. For our probe brane 
approximation to be valid, we require a large enough background charge, N > 10^. It is possible 
that if several annihilation events were to be superimposed, the amplitude of the oscillations would 
be reduced and the hierarchy between p/,/ = 1,2... become milder. Thus, we are naturally lead 
to the framework of staggered inflation, where many annihilation events are assumed to occur in 
any given Hubble time. The inclusion of many annihilation events renders a numerical study a 
necessity, which we plan to come back to in a forthcoming publication. 

To summarize, the presence of several dynamical degrees of freedom can lead to additional 
large contributions to the power-spectrum, such as a ringing pattern, if they start to decay during 
inflation. 
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